Decay of DNLS breathers through inelastic multiphonon scattering 
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We consider the long-time evolution of weakly perturbed discrete nonlinear Schrodinger breathers. 
While breather growth can occur through nonlinear interaction with one single initial linear mode, 
breather decay is found to require excitation of at least two independent modes. All growth and 
decay processes of lowest order are found to disappear for breathers larger than a threshold value. 
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The study of intrinsically localized modes in anhar- 
monic lattices, discrete breathers, has yielded much at- 
tention during the last decade (see e.g. In par- 
ticular, their existence as time-periodic solutions of non- 
linear lattice-equations was proven under quite general 
conditions |^ , and numerical schemes were developed for 
their explicit calculation [Q. The generality of the con- 
cept of discrete breathers, which provide very efficient 
means to localize energy, has lead to numerous sugges- 
tions to its application in contexts where anharmonicity 
and discreteness are important, e.g. for describing energy 
and charge transport and storage in biological macro- 
molecules [|] . Recently, discrete breathers have been ex- 
perimentally observed in coupled optical waveguides , 
in charge-density wave systems Q], in magnetic systems 
lH, and in arrays of coupled Josephson junctions [^. 

Although a discrete breather under quite general con- 
ditions is linearly stable 10 1 (and thus no pertur- 
bations grow exponentially), there are many questions 
remaining concerning the long-time fate of perturbed 
breathers. In a previous paper 0, some of these ques- 
tions were addressed considering a particular model, the 
discrete nonlinear Schrodinger (DNLS) equation. The in- 
teraction between stationary breathers and small pertur- 
bations corresponding to time-periodic cigcnsolutions to 
the linearized equations of motion around the breather 
was investigated using a multiscale perturbational ap- 
proach, and it was found that the nonlinear interaction 
between the breather and single-mode small-amplitude 
perturbations could lead to breather growth through gen- 
eration of radiating higher harmonics, but not to breather 
decay. It is the purpose of this Report to extend these 
results to more general perturbations of stationary DNLS 
breathers. In particular, we find that while the simplest 
growth process can be described as an inelastic scatter- 
ing process with a one-frequency incoming wave and an 
additional outgoing higher-harmonic wave, the descrip- 
tion of a decay process requires at least two incoming 
modes yielding outgoing modes with frequencies being 
linear combinations of the original ones. 

In order to make this Report self-contained, we first 
recapitulate the main formalism fro m ||Tl[| (to which the 
reader is referred for details; see also [|12[ for a similar ap- 
proach in continuous NLS models). With canonical con- 
jugated variables {i^/'n}, {V'nli the DNLS equation can 
be derived from the Hamiltonian 
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where A is the discrete Laplacian, (A-0)„ 
■(/)„_ 1 — 2'0„, and we will assume C > without loss of 
generality. In addition to the Hamiltonian (^, a second 
conserved quantity, the excitation norm, is defined as 



(3) 



The conservation laws for these two quantities can be ex- 
pressed in terms of continuity equations, with flux densi- 
ties for the Hamiltonian and norm given by, respectively, 



J« = -2CRe[V„+i(C+i-0 

Ja/- = 2CIm [V'rtV'ri+l 
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The single-site breather is a localized stationary solu- 
tion to Eq. (H) of the form ipn{t) = (/!)n(A)e'^*, where 
is time-independent with a single maximum at a lattice 
site (see e.g. ||l^,|,|4l). It exists for all A/C > 0, and is 
a ground state solution minimizing the Hamiltonian (|l|) 
for a fixed value of the norm (||) (see e.g. [|l^). Denoting 
the values of these quantities for the breather as Titp{h) 
(< 0) and (> 0), respectively, we have: 
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proving linear stability of the breather jl^ as well as Lya- 
punov stability for norm-conserving perturbations |14[| . 

To describe the dynamics close to the breather, we in- 
troduce the perturbation expansion (cf. ITT 



(7) 



Mt) = {0n + Ae„(t) + X''vn{t) + X'^^nit) 

-KA40„(t)-f ...}e'/^'*'. 



where e„(0) is the initial small-amplitude perturbation. 
Substituting into Eq. (^) and identifying coefficients for 
consecutive powers of the small parameter A yields for 
X^,X\\^,X^, etc: 
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etc., where (£(^^e)„ 
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Ae„. The 0th order equation (||) yields the breather 
shape {(fin} (which for the single-site breather can be 
assumed real and positive without loss of generality), 
the 1st order equation (||) is the linearization of the 
DNLS equation for small breather perturbations, while 
the higher order equations describe the long-time dynam- 
ics on increasingly longer time-scales. 

The linearized solutions are conveniently ob- 
tained by substituting e„(<) = ia ([/„ + ty„) e~*"* -t- 
^a* {U* - W*) e^'^^ into Eq. (§). This yields 
(for real (/>„) an eigenvalue problem of the form 

[C, O) [{Wn}) [{Wn}t 

operators Co and Ci defined by (cf. p^) 
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As the breather is linearly stable for all A |i5| , the eigen- 
values to are real, and the eigenvectors ({?/«}, {W^n}) 
can be chosen real and normalized. The continuous 
(phonon) spectrum of extended eigensolutions is ob- 
tained for |n| oo {(f>n 0) yielding two uncoupled 
equations for the linear combinations Un^Wn- Consider- 
ing without loss of generality only the positive-frequency 
solutions Un + Wn ^ e"^" yields the dispersion relation 
uj — A-|-2C(1— cosg), so that the continuous spectrum for 

> is the interval we [A, A -I- 4C]. Isolated eigenval- 
ues UJ ^ outside the continuous spectrum give localized 
eigensolutions corresponding to breather internal modes 
one spatially symmetric, 'breathing', mode ex- 
ists for < A/C < 1.7, and one antisymmetric, 'transla- 
tional' or 'pinning' mode for < A/C 1.1. (The varia- 
tion of their frequencies with A/C was shown in Fig. 1 of 
[ pT| ) . Finally, there are also the zero-frequency solutions 
to Eq. (^), obtained from the ansatz e„ = Un+iWn- They 
can be written as a superposition of two fundamental 
modes: the 'phase mode' |l^] Wn = 4'n describing a rota- 
tion of the overall phase of the breather, and the 'growth 
mode' |l|] [/„ = d4>n/dK (yielding e„ = d(pn/dA + i0„t) 
describing a change of breather frequency. 

This set of eigenvectors (including the zero-frequency 
modes) forms a basis for the space of solutions to Eq. (||), 
in which an arbitrary initial perturbation can be ex- 
panded. Eigenvectors (^{Un^}, {Wn^}j with different 
(real) eigenvalues w^*^ fulfill the orthogonality relations 

(w« - ^(■'")) ^ (c/Ww^i^r + = 0, (14) 



and, with the 'pseudoscalar' product defined from (|lj), 
the only nonzero product involving the zero-frequency 
modes is their cross-product p5||: 
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Let us now discuss the long-time effects of a small ini- 
tial breather perturbation, expanded in the above basis. 
The second-order correction to the linearized dynamics is 
given by the inhomogeneous Eq. (p^). As its right-hand 
side is quadratic in e, it describes fundamental scattering 
processes involving at most two modes: assuming the ini- 
tial perturbation e„(0) to consist of a set of modes with 
frequencies {w^*'}, the right-hand side will contain the 
frequencies {2wW}, {w^ -w^^']} and 0. 

Thus, it acts as a multi-periodic force localized at the 
breather region (since all terms are multiplied by </>„), 
and resonances with solutions to the homogeneous Eq. 
(§) typically yields a response for {r]n} which either di- 
verges in time (resonance with discrete spectrum) or is 
spatially unbounded (resonance with continuous part). 

In a standard way, the divergent parts can be re- 
moved by allowing for a slow, adiabatic time-evolution 
of the breather parameters. Similarly as in ||ll|, where 
the case of single-mode initial perturbations was ana- 
lyzed in detail, the divergence due to overlap between 
the static part of the right-hand side of (^o|) and the zero- 
frequency modes can be shown to be equivalent to a time- 
independent shift of the breather frequency A. Writing 
A = Ao + A^A2, where Aq is the unperturbed breather 
frequency, the second-order shift A2 will be the sum of 
the shifts resulting from the individual modes contained 
in e„(0). The latter was calculated in ||ll[ (cf. Eq. (31)), 
and generalizes for the multimode case to 

dA n i 



where a(*) is the initial amplitude for the mode 
(^{Un ^}, {Wn ^}j with frequency w*^*^. In general, A2 can 
be either positive or negative, depending on the detailed 
spatial structure of the excited eigenmodes. 

The spatially unbounded response resulting from res- 
onances between the oscillating part of the right-hand 
side of (10) and phonon modes physically corresponds to 
emission of radiation from the breather region for each os- 
cillation frequency belonging to the continuous spectrum 
of (^). The strength of the radiation fields can be cal- 
culated as illustrated in ||ll| for second-harmonic (2u;(')) 
resonances (cf. Eq. (34) and the following discussion); 
here we give the corresponding results for two-mode 
(w^'-* ± w^-'-' = w'-^-') resonances (assuming w*^*^ > w*--'-' 
without loss of generality). Writing the correspond- 



ing responses as rj- 
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are determined by 



z'^'^a'^-'-' , the functions u^^-* and li;!^-* 




i-.x r,-, u\ ■ Wiien uj'^^' belongs to 



the phonon band, the radiation field strength will be 
proportional to the overlap between this right-hand-side 
and the corresponding continuous spectrum eigenvector 

({Uk^^}, {Wi^'}) , which using (HI) is obtained as 



n 



(16) 



Far away from the breather, the radiation field should 
correspond to two identical outgoing propagating linear 
waves, yielding the boundary conditions wi^'' — > 



arccos 1 



(i-»^) 



±00 where r^^^ 



c(±) and = 

from the linear dispersion relation. 

Let us now discuss the consequences of the second- 
order radiation arising from the two-mode- interaction for 
the breather itself. Assume for simplicity that only two 
modes with frequencies w'-"'^-' and w^^' {uj'^^^ > w^^^) are 
initially excited (the presence of other modes will con- 
tribute to orders 3 and higher), and that both modes 
belong to the continuous spectrum (the case of internal 
mode excitation requires some modification, see p| . Sec. 
IIIC). Far away from the breather, we assume a station- 
ary regime which, in the most general case when 2iLj^^\ 
2(jj'^', Lj^^^ and a;'-^-' all belong to the phonon spectrum, 
corresponds to the boundary conditions at n — > ±c»: 
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+a(2)(eT^«'''" + r(2)e±^«'''") 

+ (a(i))2r«e^(±9^''"-2-'^'*) 

-fa(i)(a(2))v(-)e'(±«'"'"-"'"'*) 



(17) 



Note that in general, the stationary amplitudes for the 
outgoing fundamental waves will differ from those of the 
incoming (i.e., r*^^^ , r*^^^ 7^ 1); this is a consequence of res- 
onances at the frequencies w'^' and w'^' in Eq. (^Tj) for 
the third-order field ^„ . For a general multimode pertur- 
bation, it is seen from (^W that this correction takes the 
form r^*' = 1 -|- X)j 1^"''' where the sum goes over 
all initially excited modes (here j = 1, 2). 

We then consider the conservation laws for the total 
norm and Hamiltonian, respectively, contained in some 
region around the breather. Assuming the breather fre- 
quency A to be the only time-dependent parameter in 



the stationary regime we can, similarly as in Ref. 
Sec. IV, write the time-averaged balance equations as 
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dt 
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A = {J^fi~^))t'{J^i+'x))t, (18) 



A= (J„(-«)))t-(Jw(+oo))t. (19) 



As the time-average of the flux densities J/v and J-h 
are additive quantities for small-amplitude plane waves 
of the form ipn — yle**^'^""^*^ the right-hand sides are 
readily obtained from ( p^ using the general expressions 
Jj^ = 2|ApCsinQ and J-h — VlJj^ for the individual 
waves resulting from (|)-(|) Combining Eqs. (|l|)-(|l|) 
and using (^ yields the time-derivative of the breather 
frequency to order 4 in the mode amplitudes a(i),a(2) as: 
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+ |a(2)p(l-5|)(l-p)p)sin,(2) 
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For a single-mode initial excitation, a^^-* = and only 
the first, positive term of the right-hand side remains, 
so that we recover the result (Eq. (55) in ||ll|) that the 
generation of second-harmonic radiation always leads to 
breather nrowth. However, in the two-mode case, the ex- 
pression (20) is in general not sign-definite, and in partic- 
ular the contribution from the frequency w^^-* is always 
negative. It is therefore natural to associate this radi- 
ation with the fundamental lowest- order mechanism for 
breather decay. This can be seen as a consequence of 
the Hamiltonian flux density being proportional to the 
frequency for plane waves. Thus the scattering towards 
the lower frequency w'"^ should yield a net flow of neg- 
ative Hamiltonian energy out from the breather region, 
to which the breather adapts by decreasing its frequency 
and amplitude according to (|6|). On the other hand, the 
generation of radiation with the higher frequency cj^"*"^ 
should analogously contribute to breather growth, as for 
the second-harmonic case. (See also for a similar 
explanation of breather decay in a Klein-Gordon model 
resulting from a resonance in the linearized equations.) 

<"a + 4C, 

so that we have < w'^' - w^^) < 4C and 2A < 
2u;(i\2a;(2),a;(i)+w(2) < 2A-h8C. Thus, we see that sec- 
ond order radiation can only be generated if the breather 
frequency fulfils < A < 4C, since for A > 4C the 
frequency tj*^^^ is always below the phonon band, while 
2w(i),2tj(2) and are always above. Therefore, when 



It is important to note, that when w^^-* and oj^"' are 



phonon modes, they fulfill A 



3 



A > AC all terms in the right-hand side of Eq. (|2C 
necessarily vanish, and all growth and decay processes 
of lowest order disappear. Numerically, we found g 
that this corresponds to breathers with central-site in- 
tensity IV'noP ^ 5.65. Therefore, these large-amplitude 
breathers are particularly stable, since all possible growth 
and decay processes result from third- and higher-order 
radiation processes so that the rate of their growth/decay 
must be at least of order 6 in the initial mode amplitudes. 

To illustrate the fundamental lowest-order decay mech- 
anism, we consider the 'pure' case when uj^~'> belongs to 
the phonon spectrum, while 2uj'^^\ 2uj'^-^\ and ijJ*^+^ are 
outside. Thus, we have ro^-* = ro^-* = r^^'> ~ 0, and 



only the last two terms in Eq. ( |20|) are nonzero. We note 



that if jr^^'l < 1, the last term will be positive, and we 
can therefore in general not conclude that A must be 
negative. The reason is, that although the scattering to- 
wards the lower frequency w^"-* from either of the two 
frequencies w*-^^ or w^^-* considered independently would 
yield a net outflow of negative Hamiltonian energy from 
the breather region and thereby breather decay, there 
will also be a contribution from the mixing between the 
Lo^^^ and oj^^^ modes arising from the third order equa- 
tion This contribution would yield a net outflow of 



positive Hamiltonian energy if 



r(2) 



< 1 and |r(i)| > 1, 



and could therefore contribute to breather growth. How- 
ever, extensive numerical simulations for different A, o;*-^^ 
and w*-^^ belonging to this region of 'pure' ^^"^-scattering 
have always shown that the net result is a linear decrease 
of the breather frequency, and therefore we believe it jus- 
tified to identify this two-wave scattering as the funda- 
mental lowest-order breather decay mechanism (Fig. |l|). 
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FIG. 1. Time evolution of the total norm M contained in a 
region of 120 sites around a breather with frequency A = 0.45, 
perturbed along two eigenmodes with frequencies o;^^' ~ 3.00 
and tj'^' 2.47, respectively {C 

In conclusion, we found that while in the DNLS model 
breather growth can result through interaction between 
the breather and single-mode initial perturbations, the 
description of breather decay requires simultaneous ex- 
citation of at least two independent linear modes. This 



confirms numerical results in Ref. (Figs. 5 and 6), 
showing breather decay also from initially single-mode 
perturbations with larger amplitude, as more frequencies 
became gradually excited e.g. through oscillatory wave 
instabilities We believe that our approach could be 
useful to understand the properties of the stationary in- 
tensity probability distribution function in the 'negative 
temperature'-regime of the DNLS model [Q, where per- 
sistent localized breathers were found, weakly interacting 
with small amplitude radiation. Finally, we remark that 
the DNLS model is nongeneric among Hamiltonian lat- 
tice models as it has two conserved quantities, and the 
DNLS breather has only one fundamental frequency with 
no harmonics. Thus, our approach cannot immediately 
be extended to other models exhibiting breathers such as 
Fermi-Pasta-Ulam or Klein-Gordon lattices. However, 
as for the latter the DNLS equation is known (see e.g. 

to describe the small-amplitude dynamics for small 
inter-site coupling, we believe that the breather growth 
and decay mechanisms described here are relevant also in 
these systems. This will be investigated in future work. 
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